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Abstract 

A scenario for laser-driven ion acceleration from a thin solid target is considered, 
where the reflected laser pulse is redirected towards the target by reflection at the 
inner cavity wall. This scenario is discussed in the context of sub-wavelength 
foil acceleration in the radiation pressure regime, when plasma dynamics is 
known to be reasonably well described by the laser-sail model. A semi-analytic 
extension of the ID laser-sail model is constructed for this case. The effect 
of cavity reflections on sub-wavelength foil acceleration is then illustrated with 
two concrete examples of intense laser pulses of picosecond and femtosecond 
duration. 
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1. Introduction 

In the last decade there has been growing interest in laser-driven ion ac- 
celeration, where energetic ions are generated by irradiating a thin foil (with 
thickness in the micrometer range) with an ultra-intense (above 10 18 W/cm 2 ) 
laser beam. The laser-generated proton pulses had been used for example for 
ultra-fast radiography and radioisotope production, and together with laser- 
generated beams of heavier ions they are of high interest from the point of 
hadron tumor therapy and fast ignition of inertial confinement fusion targets 
[H, 0, 0, 0] ■ The use of an ultra- intense laser as an ion driver offers certain advan- 
tages over conventional linear or circular accelerators, allowing for construction 
of very compact devices that produce particle pulses of high intensity and very 
short duration. Unfortunately, ion energies achieved so far in this approach 
are not very high (below 100 MeV per nucleon), and the achieved values of the 
laser-to-ion energy conversion efficiency are generally below 10% [H, 0|. 
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Figure 1: Schematic drawing of a target with the cavity that redirects the reflected pulse 
towards the foil. 



It was recently observed that basic parameters of the accelerated ion pulses 
- the maximum ion energy, the average ion energy, and the laser-to-ion energy 
conversion efficiency - could be considerably improved if the target foil is en- 
closed in a cavity that redirects the laser pulse reflected from the foil back into 
the foil, thus in a sense "recycling" the laser energy [7, 8]. The mechanism of 
cavity- enhanced plasma acceleration, dubbed the Laser Induced Cavity Pres- 
sure Acceleration (LICPA) mechanism H Q , was found to be very effective at 
sub-relativistic laser intensities (~ 10 15 W/cm 2 ), where enhancement in ion ac- 
celeration ispredominantly due to the buildup of thermal plasma pressure inside 
the cavity 0, @, E3] ■ Computer simulations performed at higher laser intensi- 
ties, when the radiation pressure becomes important , show that the LICPA 
mechanism would be effective in this regime as well 0, HI]. A possible target 
arrangement realizing the LICPA principle in the radiation pressure regime is 
shown in Fig. [TJ A detailed account of these simulations would be reported 
elsewhere. In this note we discuss a realization of the LICPA idea in the radia- 
tion pressure regime, but in a narrower context, namely when plasma dynamics 
is reasonably well described by the laser-sail (LS) model (ill flil fl3 | . Such 
a situation occurs for example when a high intensity and high contrast laser 
pulse is incident on a foil of sub- wavelength thickness. This acceleration regime 
attracted recently considerable attention [3 [H, [H, [l7|, El, [H, S ED]. We 
limit our considerations to a one-dimensional approximation and we construct 
cavity- enhanced LS model, in which the effect of reflections of the laser pulse at 
the inner cavity wall is taken into account in a semi-analytic way. We illustrate 
the characteristic features of this implementation of LICPA with concrete pre- 
dictions for ions generated using a realistic ultra-thin foil target and the laser 
parameters close to the expected capabilities of the big new laser facilities ELI 



2 



22 1 and PETAL [23| that would become available in the near future. 



2. The laser-sail model in one dimension 

Let us consider a thin planar foil of areal mass density a, irradiated at a 
normal incidence by an intense laser beam propagating along the x axis. The 
intensity of the laser beam at the origin of the laboratory reference frame has 
the temporal profile I(t). Under certain circumstances - for example when 
the foil has sub-wavelength thickness and the laser pulse is of sufficiently high 
intensity and high contrast - the overall motion of the accelerated plasma is 
similar to that of a "sail" driven by the radiation pressure of the laser beam 
(l3j ). Furthermore, if the foil thickness exceeds certain minimal value, it acts 
as a nearly totally reflecting mirror. The equation of motion for such a mirror 
driven by the laser beam had been first derived in the context of laser-driven 
interstellar space travel [HE3: 



« m = wt^M RW) ^ (1) 

at oc l 1 + p 



where P = 7 = l/^l — P 2 , lu' = — PT/O- + P) is the laser frequency 

in the foil rest frame, and R(u)') denotes the foil reflectivity in its rest frame. We 
shall further assume R(lo') = 1. This equation may be recast in a particularly 
useful form by introducing a relativistic retarded time variable w = ct — x fl2| : 

7 dp _ 2 I(w/c) _ (2) 



1 — p dw (TC 3 

By a straightforward integration it is then found that the speed of the foil which 
was at rest at w — is given by 

(l + e( W )) 2 -f 

where e(w) is the total laser energy per unit area incident on the foil up to the 
retarded time w, normalized to the half of the relativistic rest energy per unit 
area of the foil: 

2 f w 

e(w) = / I(w'/c)dw' . (4) 
ac Jo 

Using the fact that dx/dw = P/(l — P) one may obtain the position of the foil 
as a function of the retarded time w: 

x(w) = x(0) + J™ Y^iu/) dW> ' ^ 

Taking into account that ct(w) = w + x(w) we obtain this way a complete 
parametric description for the world line of the foil in the (ct, x) plane in a 
parametric form. The relativistic kinetic energy of the ions in the foil is found 
to be 
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This very simple ID model was found to be surprisingly useful in interpret- 
ing results of 3D and 2D kinetic simulations for laser ion acceleration. In the 
case of sub-wavelength foils driven by linearly polarized beams, the laser sail 
behaviour was famously identified in a 3D simulation at the laser intensity of 
1.4 x 10 23 W/cm 2 It was found that the initially planar foil is deformed 

into a "cocoon", in which the laser pulse is confined, and a stable plasma clump 
is formed from the portion of the foil of the size of the pulse focal spot, being 
then pushed forward by the radiation pressure. (In fact the LS model was used 
in that paper to extrapolate the ion acceleration beyond the simulation time, 
which was limited due to insufficient computer resources.) It was later observed 
that circularly polarized beams create more favorable conditions for such a block 
acceleration of plasma [24], and a parametric study had shown that the LS-type 
of behavior becomes visible in ID simulations for ultra-thin targets already at 
I ■ (A/ljUm) 2 = 2.8 x 10 20 W/cm 2 [2(J. The case of linear polarization is more 
complex, because generally at less extreme intensities a broad energy spectrum 
is obtained. The most we can hope for is a quasi-monoenergetic peak on a broad 
background, which becomes clearly visible at the intensity of 5xl0 21 W/cm 2 



25l . l2Gl | . It seems then reasonable to expect that evolution of ions in the quasi- 
monoenergetic peak would be consistent with the LS model. It is interesting to 
note that the evolution of plasma clumps follows the LS dynamics even though 
they are not accelerated as rigid bodies and undergo some internal evolution. 
It is also interesting that the evolution of targets which are relatively thick 
compared to the laser wavelength is also found to be consistent with the LS dy- 
namics (27, 28 1 . A gen eralization of the relativistic LS foil dynamics beyond ID 
was outlined in fbH Il8| and applied to the study of instabilities in the ultra-thin 
foil acceleration process. 



3. Cavity-enhanced laser-sail model 

Let us now consider a generalization of the ID laser-sail equation to include 
a contribution from the laser pulse reflected from the target and redirected onto 
the target through a reflection from the inner wall of the cavity surrounding the 
target. The enhanced accelerating effect of the laser pulse reflected from the 
target foil and then redirected towards the target was considered already in fl2l | 
in the context of an interstellar space travel driven by an Earth-based laser, but 
for the laser ion acceleration a more general formula is needed, valid also in the 
case of an overlap between the original pulse and the pulse reflected at the inner 
cavity wall. 

Let us assume that the laser pulse enters the cavity at x — and is totally 
reflected by the foil located at x(w) (this implies that the foil thickness - being 
sub- wavelength - has to be larger than certain minimal value) . The pulse returns 
to the inner cavity wall located at x — and is then redirected towards the foil 
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(this is of course an approximation, in a real 3D geometry the reflection at the 
inner cavity wall would have to occur in a different place than the entrance to 
the cavity). The distance between the cavity entrance at x = and the initial 
position of the foil x(0) is equal to the cavity depth L c . Let us denote by Wj, 
j = 0, 1, 2...., the values of the retarded time for which the leading edge of the 
laser pulse is sent toward the foil from the position of the inner cavity wall for 
the j-th time. Assuming that the laser pulse begins at wq = 0, we obviously 
have W\ = 2L C . Next, let x^(w) denote (for j > 1) the position of the foil 
for w in the interval [wj—±,Wj]. These functions satisfy an obvious continuity 
relation x^ +1 ^ (wj) = x^\wj). Knowing Wj and x^(wj) , we obtain Wj+x from 
the formula 

w j+ i = Wj +2xW(wj) . (7) 

Finally, let e^\w) denote (for j > 1) the total laser energy per unit area of 
the foil - expressed in the units of a half of the relativistic rest energy per unit 
area of the foil - that was incident on the foil for w in the interval [wo,Wj], i.e. 
during the whole acceleration process up to Wj. The value of e^(w) is simply 
given by the formula 2J The values of e^(w) for j > 2 may be determined 
iteratively. The expression for e^ +1 \w) consists of three terms: the first repre- 
senting the total laser energy per unit area incident on the foil until w — Wj, the 
second representing the laser energy per unit area coming from the original laser 
pulse, incident on the foil in the interval [wj , w] , and the third representing the 
laser energy per unit area that had been reflected from the foil in the interval 
[tOj—i, Wj] and was redirected towards the foil. The laser energy per unit area 
reflected from the foil in an interval [w a , Wb] is given by the formula 

ac 2 e(w b ) - e(wg) 
2 [l + eK)][l + e«)] ' 1 ' 

but only a fraction R c of this is redirected to the foil by the reflection at the inner 
cavity wall. (For simplicity we shall assume in the following that R c is constant, 
whereas in reality it could show some dependence on laser intensity and material 
properties of the inner cavity wall.) Taking into account all contributions we 
obtain 



U+Dt ^ U)f ^ 2 r u '/w D e^{w% {w))-e^{wj. x ) 
^ +l '(w) = e w (wj)-\ — / I(w /c)dw +R r 



[l + eU)( w ^ f (w))][l+eU)( Wj _ 1 )] , 

(9) 

where the function w^t {w) is a solution of an implicit equation 

w = w% + 2x^{w%) 1 (10) 

i.e. it represents the value of the retarded time variable from the interval 
[v)j—i, Wj] for which a redirected ray that is incident on the foil at w £ [Wj , tfty+j.] 
was for the last time reflected from the foil. Obviously, w^Awj) — Wj-\. The 



relativistic speed f3^ +1 \w) in the interval w 6 [wj,Wj+i] is obtained by substi- 
tuting e^ +1 ^(w) into the Eq. [3l The position of the foil is then obtained from 
the formula similar to Eq. [SJ 

The relativistic kinetic energy of the accelerated ions for w € [wo,Wj] is given 
by the formula |6] with e(w) replaced by e^'(w). 

4. Concrete examples 

Let us consider the predictions of the cavity-enhanced LS model in two con- 
crete scenarios characterized by the same total laser fluence, but with distinctly 
different pulse lengths. The first scenario involves a laser pulse with the max- 
imum intensity of 10 22 W/cm 2 and the full width at half maximum equal to 
tnw = 130 fs; such parameters should be attainable for example with a slightly 



defocused kilojoule beam at the future ELI Beamlines facility 22]. The sec- 
ond scenario involves a pulse with the duration of 1.3 ps, i.e. larger by an 
order magnitude, but with lower maximum intensity, equal to 10 21 W/cm 2 ; 
such parameters should be available for example at the future PETAL facility 
[23| | . The temporal profiles of these pulses are taken to be super-Gaussian, and 
they are assumed to start at the intensity smaller by a factor of 100 relative 
to the maximum intensity: I(t) — I exp[—((t — < p /2)/r) 6 ] for t € [0,t p ], where 
t p = twQoglOO/log^) 1 / 6 and r = ±t HW /(log2) 1 / 6 . The laser energy fluence 
on target in both scenarios is equal to 12.81 J/(/im) 2 . 

Our target is assumed to be an ultra-thin carbon foil, with the density of 
2.2 g/cm 3 . For the foil to be nearly completely reflective it cannot be too 
thin. In ID the condition for nearly total reflection takes the form (20L l2l| 
7r(7/A) • (n e /n c ) > ao, where A is the laser wavelength (which is either 800 nm 
for ELI or 1054 nm for an Nd:glass picosecond laser), n e is the plasma electron 
density, n c = (4ir 2 c 2 m e Eo/e 2 ) ■ (1/A 2 ) is the critical electron density for a given 
wavelength (m e is the electron mass), and ao is the dimensionless laser field 
amplitude, related to the intensity / by the relation / = otm e (?ri c a^ {a. = 1 for 
circular polarization and a — 1/2 for linear polarization). The non-transparency 
condition may be rewritten in the form 

(12) 

For a fully ionized carbon we have n e — 6.62 • 10 23 cm" 3 , so for a circularly 
polarized beam with / = 10 22 W/cm 2 we obtain a condition I > 32 nm. How- 
ever, to leave some margin for the foil thinning due to the 3D stretching, as well 
as possible prepulse issues which might be important in the case of kilojoule 
beams, we assume / = 200 nm. 

As for the cavity parameters, in both cases we choose L c = 80 /im, which 
seems to be a technologically achievable value. The reflection coefficient at the 
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Figure 2: Ion energy per atomic mass unit, as a function of the acceleration distance, for 
the case of 1.3 ps pulse with 10 21 W/cm 2 peak intensity. Solid line indicates acceleration 
without cavity reflection, while the dash-dotted (dashed) line indicates an acceleration with 
the fraction R c = 0.8 (R c = 0.6) of the laser energy being reflected at the inner cavity wall. 

inner cavity wall R c is a free parameter in this model. We assume R c = 0.8 and 
for comparison we give also the results obtained for R c = 0.6; a more precise 
determination of R c would require computer modelling of the cavity beyond the 
ID approximation. 

In Fig. [2] we show the evolution of the ion energy per atomic mass unit 
of ion mass (i.e. E^/A, where A is the ion mass number) as a function of 
the acceleration distance for the 1.3 ps pulse with 10 21 W/cm 2 peak intensity. 
The evolution in the absence of the cavity reflections is indicated by the solid 
line: the foil is rapidly accelerated over a distance of 215 /jm to the energy of 
119 MeV/u, and then "sails away" at a constant speed. 

If reflections at the inner cavity wall are taken into account with the fraction 
of R c = 0.8 of the laser energy reflected at the inner cavity wall , the acceleration 
of the foil is continuously enhanced by the laser pulse redirected multiple time 
towards the target by reflection at the inner cavity wall (as indicated by the 
dash-dotted line), which nearly doubles the resulting ion energy relative to the 
case without the cavity reflection. For R c = 0.6 (dashed line) the energies to 
which the ions are accelerated are slightly smaller, but the general pattern is 
the same. 
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Figure 3: A space-time plot of successive reflections of the incident laser pulse from the 
accelerated foil and the inner cavity wall for the 1.3 ps pulse with 10 21 W/cm 2 peak intensity, 
with the fraction R c = 0.8 of the laser energy being reflected at the inner cavity wall. 



The result for the ion energy per a.m.u. is easily converted into the laser- 
to-ion conversion efficiency, which is shown on the right side of the figure. The 
total ion energy per unit surface area of the target is given by {E\™jA) ■ (cr/u), 
where u is the atomic mass unit, so every 10 MeV of ion energy per a.m.u. 
corresponds to the ion beam energy fluence of 42.45 MJ/cm 2 , which implies a 
laser-to-ion energy conversion efficiency of 3.31%. Note that the acceleration 
efficiency is quite high - around 40% - even without the cavity reflections, which 
is characteristic of the radiation pressure acceleration 2(| • 

A space-time plot of successive reflections of the incident pulse in this case is 
presented in Fig. [3] The thick segment along the vertical axis indicates duration 
of the initial pulse at the entrance to the cavity, the moderately thick solid curve 
indicates the world line of the foil, the thin solid line indicates the world line of 
the leading edge of the pulse, including the reflections from the accelerated foil 
and the inner cavity wall, and the dashed line indicates the world line of the 
trailing edge of the pulse. As we see, there is a substantial overlap between the 
incident pulse and the reflected pulse in this case. 

The evolution of ion energy per a.m.u. for the case of 130 fs pulse with 
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Figure 4: Ion energy per atomic mass unit, as a function of the acceleration distance, for 
the 130 fs pulse with 10 22 W/cm 2 peak intensity. Solid line indicates acceleration without 
cavity reflections, while the dash-dotted (dashed) line indicates the acceleration with the 
fractionR c = 0.8 (R c = 0.6) of the laser energy being reflected at the inner cavity wall. 



10 22 W/cm 2 peak intensity is shown in Fig. |U If there are no cavity reflections, 
the foil is rapidly accelerated over a very short distance of 22 /im, achieving the 
same energy as in the previous case, and then "sails away" at a constant speed. 
If reflections at the inner cavity wall are taken into account (with R c = 0.8), 
the foil receives additional kick by the redirected pulse, which re-interacts with 
the foil at the distance of 113 from its initial position, again nearly doubling 
the ion energy (and the conversion efficiency as well). After being reflected by 
the foil for the second time the pulse is redirected onto the foil once again and 
hits the moving foil for the third time at the distance of 645 /im from its initial 
position, but the effect of this third impact is rather small. It should be noted 
that the spatial extent of the femtosecond pulse is small compared to the chosen 
cavity depth, so the foil flies freely before being hit by the redirected pulse. For 
R c = 0.6 the energies of the accelerated ions are slightly smaller, but the general 
pattern is the same. 

Interestingly, energies achieved over the acceleration distance of several hun- 
dred micrometers in this scenario are very close to energies achieved with the 
ps pulse, which indicates that higher order reflections have only small effect 
on the resulting ion energies. This is an important observation from the point 
of view of designing targets with a cavity, because it shows that it is enough 
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for the cavity reflector to survive only few reflections, since the loss of further 
reflections due to a deteriorating cavity wall would not significantly affect the 
final result. 

However, this picture could be different if a more dense target is used, for 
example a metal foil. Such foil would accelerate less rapidly, so even for a small 
acceleration distance many cavity reflections could occur. The ion energy per 
a.m.u. achieved in the pure LS mode would be small in this case, but the 
enhancement due to presence of the cavity could be much bigger than 100%. 
The question is whether the foil and the cavity would survive as a solid body 
after large number of reflections. 

Finally, it is an interesting issue how the acceleration energies depend on 
the depth of the cavity L c . We studied this issue for the case of the picosecond 
pulse. Somewhat surprisingly, it was found that the ion energies measured at 
the acceleration distance of 300 ^m show only a mild dependence on the cavity 
depth: for L c = 40, 60, 80, 100, 120, 140 /im the values of ion kinetic energy per 
a.m.u. were respectively 208, 202, 195, 186, 179, 175 MeV. For the femtosecond 
pulse the dependence would be even smaller, since even for the smallest of the 
considered cavities there would be no overlap between the incident pulse and 
the redirected pulse. 

5. Conclusions 

Summarizing, we constructed a simple implementation of the LICPA mech- 
anism for laser ion acceleration in the radiation pressure domain, based on a 
semi-analytic ID laser-sail model. The model accomodates in a general way the 
additional contribution from pulses reflected from the target foil and redirected 
towards the target via reflection at the inner cavity wall. We used this model to 
obtain concrete predictions for acceleration of an ultra-thin carbon foil, driven 
by laser pulses of different duration - a f s range pulse and a ps range pulse - and 
different intensity, but the same laser fluence, providing a convincing illustration 
of characteristic features of LICPA. The possibility of a considerable increase 
in the ion energy and the laser-to-ion conversion efficiency, as large as 50% or 
more, was clearly demonstrated. Our predictions could serve as a guide for ex- 
periments at the laser facilities that would become available in the near future. 
Although in principle the cavity-enhanced LS model is applicable to ultra-thin 
foil acceleration, it is not unreasonable to expect that it would turn out to be 
useful in interpreting the results for averaged quantities obtained in numerical 
simulations and experiments involving thicker targets, as had been indicated in 

P. R^czka gratefully acknowledges ESF support in the form of the SILMI 
RNP grant No. 3638. 
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